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I. INTRODUCTION 



Interactions of the rj meson with few-nucleon systems complement our knowledge of the 77- 
nucleon interaction. An independent, although related, interest in these systems stems from 
the conjectured existence of 77-nuclear quasi-bound states. Such states have been predicted 
by Haider and Liu |IJ, and Li et al. @, when it was realised that the 7/-nucleon interaction 
is attractive. So far there has been no direct experimental verification of this hypothesis 
but, on the other hand, there is a mounting evidence that such states exist in the lightest 
few-nucleon systems. 

The earliest measurement was done in the four-body case. It was found at SATURNE 
that the pd — > r/ 3 He production amplitude falls rapidly just above the rj production threshold, 
||, j§. This led Wilkin to the suggestion that such a quasi-bound state has been produced 
0. A similar but weaker slope was found in the dd — > r/ 4 He amplitude ||, and this may 
indicate that the quasi-bound states in five-body systems are broader and located further 
away from the rj meson production threshold |J. 

In more recent studies of three-body systems, very strong final state interactions (FSI) 
were found in the pp — > pprj cross section in the threshold region [RJ, H, M. Those led to 



a speculation that this system may form a quasi-borromean state, ||10|| . As the final state 
interactions in the three-body continuum are difficult to analyse, one finds perhaps more 
convincing evidence in the measurement performed at CELSIUS, [PH| . There, the 7/-deuteron 
final states are observed in the scattering of a proton on a deuteron target. This method 
allows the determination of the pn —>■ rjd amplitude, which turns out to be a smooth function 
of energy but displays an enhancement close to the threshold. Much stronger enhancement 
is indicated by an earlier SATURNE measurements performed by neutron scattering on a 



proton target [12]. These data taken together indicate a very sharp threshold effect in the 



7/-deuteron system. Earlier calculations performed by Ueda confirm such an effect ||13|| , and 
indicate the existence of quasi-bound i]NN states. 

The correlations seen in the three-body systems may be due to quasi-bound states or 
quasi-virtual states, i.e. three-body analogs of the deuteron or the virtual spin singlet state 
in the NN system. However, there has been a serious difficulty in such an interpretation. 
The strength of the r]N interaction and the magnitude of the scattering length a V N were 
not large enough. The early values of Re a V N ranged between 0.25 fm and 0.5 fm , |T|, [Q, 



0, [|15|1 , |16| , |jj although larger values have also been suggested |18fl , |19| , pCf] , |2T 



|22[. The lengths up to 0.5 fm are not large enough to support long lived structures in 
rjNN systems. Calculations have shown that the ^-deuteron virtual state is only likely to 



be formed provided Re a V N exceeds a critical value of about 0.7 fm , [p3[ . For larger values 



of about 1.2 fm, one finds the rjd system likely to be quasi-bound p3| . Such states would 
be reflected in a large value of the rjd scattering length A v d- Unfortunately, the calculations 
of A v d done in this critical region depend on details of the r]N interaction model. The 
significance of the scattering length is one factor, but another equally important one is the 
way the r]N scattering amplitude extrapolates to the region below the r]N threshold. 

Recent phenomenological analyses of a four coupled-channel (t]N, irN, ^N, ttitN) K- 
matrix model have helped to resolve these difficulties. In Ref. ]25| a large scattering length 
0"qN = 0.75(4) + i0.27(3)/m was found, and in addition the effective range expansion was 
shown to work well in a broad region of both positive and negative energies. A very similar 
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value of a V N = 0.72(3) + z0.26(2)/m was obtained in a dispersion theoretic analysis of 
reference [f26|| . However, further studies [27] of the K-matrix model done in the spirit of 
Ref. [25[] have now revealed the existence of another solution, which is characterised by a 
small scattering length of pa 0.3 + i0.2fm and a very large effective range. It should be 
added that a low value of a is also produced by the GW21 solution in the authors iV(1535) 
pole-position studies of Ref. | 28| . [] In that work only the channels 7rN, rjN and tcttN were 



treated explicitly. However, when we use an extended version of the .ff -matrix model, in 
which both the jir and 777 channels are included |29| then this GW21 solution disappears 
leaving only GW11 with its large scattering length Re a rj ^ = 0.87/m. One hopes that the 
few-body data are going to be selective on these possibilities. In addition, if a vN is large 
there should arise a rich spectroscopy of nuclear 77 states. Forthcoming experiments may 
resolve this question, PD[ , although difficulties in the interpretation have been suggested, 

It seems that now there exist enough information to study the rjNN phenomena in an 
almost quantitative way, and some related problems are studied in this paper : 
In section II, the 77-deuteron scattering amplitude at low energies is calculated. A method 
of partial summation of the multiple scattering series used previously to calculate the eta- 
deuteron scattering length, |2"3|] , is extended to the case of the continuum. A strong cusp in 
the rjd scattering amplitude at the threshold is found and attributed to a virtual state. The 
width of this state is very small, in the 1 MeV range. 

In section III, the final state interactions in the rjd system are discussed. We find the virtual 
state to couple rather weakly to this system. It produces rather limited effects at energies 
below the deuteron breakup. The results of Ueda, fll3| , who has predicted a broad and 
prominent peak above the rjd threshold, are not reproduced. 

In section IV, we discuss the constraints induced upon the rjN scattering matrices by the 
rjNN data. Also, the nature of exotic rjNN states is discussed in this section. 

Some parts of this research could be, and have been, done in an exact way by solving the 
three-body Faddeev equations plfl , [[3^], f3"3fl . The solution for rjd scattering is the simplest 
one that can be improved in this way. However, even for this simple question the general 
multiple channel analysis has never been fully implemented. It seems also to be premature, 
as many ingredients of the rjN interaction are not known. The advantage of the method 
presented here consists in the following : 

(1) It is numerically simple and transparent in its physical interpretation. Uncertain- 
ties in the two-body interactions and the meson formation process may be separated and 
discussed in simple terms. 

(2) This method is applicable to more complicated few nucleon systems. 

(3) It is very accurate numerically, at least in the calculation of the rjd scattering length. 



The authors wish to thank Bengt Karlsson for emphasising this point. 
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II. THE 77 DEUTERON SYSTEM AT LOW ENERGIES 



We discuss the rjd scattering at very low energies. This section contains three parts. First 
a simple method is introduced which allows us to find the scattering matrix T^. At very low 
energies, below deuteron breakup, T vd is dominated by a narrow virtual three-body state. 
Next, we study the off-shell extension of T v d in order to construct the r]d wave function. The 
latter is used to describe final state interactions in the rjd formation processes. In the third 
part the static nucleon approximation is discussed. Everywhere throughout this paper the 
dominance of S waves is assumed. 



A. A Formula for the r/-deuteron T matrix 

The purpose of this subsection is to derive a simple formula which relates the low-energy 
meson-deuteron scattering amplitude to a few multiple scattering integrals involving the 



meson-nucleon amplitudes. The formula is described in more detail in Ref. [E3J, where it 



was used to calculate the rjd scattering length. Its numerical precision was later confirmed 



by a direct solution of this problem in terms of the Faddeev equation . In the limiting 
case of fixed nucleons an agreement with other calculations has also been established [38 . 
Here, this simple procedure is outlined and its extension to the elastic or inelastic continuum 
and to the off-shell region is discussed. 

The starting point is a fixed nucleon approximation. Next, upon the solution of this 
simple limit a more realistic calculation is built up. It is achieved in a step-by- step summa- 
tion of the multiple scattering series until the necessary precision is reached. The scattering 
length of a meson on a pair of fixed nucleons is given by 

_ 2a vN 

vd ~ 1 _ a v N ' ^ > 

Rd 

where a V N is the meson-nucleon scattering length and Rd is the inter-nucleon distance, [BU 



Eq.([l]) is obtained in a simple way by setting boundary conditions for the meson wave 
function if) at each scatterer as ip'/ip = l/a V N- An alternative way to obtain the same 
result is to sum the multiple scattering series, which in this case is a geometric series. With 
this expansion, the l/Rd term stands for the meson propagator while a V N j^^n is a double 
scattering amplitude etc. The main advantage of Eq.([l]) is that it offers the solution even 
for a divergent multiple scattering series when / Rd > 1- The main weakness is that it is 
based on distinguishable nucleons. Also, important kinematical and dynamical effects of the 
three body propagation are missing and the actual value of Rd is left undetermined. The 
method to be described below retains the main advantage, but removes the weak points and 
offers a correct value of Rd- 

To improve Eq. (|I]) for the non- fixed nucleon situation one needs to introduce: the three- 
body propagator, the NN interactions and an off-shell t v ^ scattering matrix. To achieve this, 
we study the multiple scattering expansion for meson-deuteron scattering. The basic NN 
interactions are assumed to be known and described by the deuteron wave function (pd and 
a complete off-shell Inn scattering matrix. With t^N given one finds a partial three-body 
propagator 
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G — Go + Gnn, (2) 

where Go is the free three-body propagator and 

Gnn = Got^NGo (3) 

constitutes that part of the three-body propagator which describes the NN interactions. 
Now, the series appropriate for the 7/-deuteron scattering matrix T vd follows from the Faddeev 
equations : 

T v d = ti +t 2 + t\G()t2 + t 2 Got\ + t 2 Got\Got 2 + tiG^Goti 

+ (ti + i 2 )G NN (t 1 + t 2 ) + (k + i 2 )G NN (k + f 2 )Gjvjv(ti + k) + 

where the nucleons are labeled 1,2 and the ti are ^-nucleon scattering matrices. This ex- 
pansion is performed in momentum space. All quantities involved in Eq.@ are operators 
and this is denoted by the hat Q symbols. Thus, appropriate multifold integrations over the 
intermediate Jacobi momenta (p, q) must be performed. In full generality the kernel of T v d 
is a four point function T vd (p, q, q',p'; E) and to find the r/d scattering matrix T vd a matrix 
element has to be calculated 

T vd (pi,E,p f ) =< <j) d ?Pri(Pi) I T vd | (f>dipr,{Pf) >, (5) 

where the ip v are S- wave functions for the initial and final mesons of momenta Pi, Pf- In 
this section we are mainly concerned with the on-shell value of this scattering matrix, which 
is denoted by T v d(E) with E being the nonrelativistic energy of the rjd system. The leading 
order (impulse approximation) matrix is 

T° d =<f° d >=<i 1 + i 2 >, (6) 

where <> denotes the average over the deuteron and mesonic wave functions as indicated 
in Eq.([|). The first partial sum of the series for T vd (E) is obtained with the expression 



< > 
1 - Si - fii 

which contains the double scattering terms 



Tid = , ^ n (7) 



< T® d G NN T® d > 

1 = <^> (8) 



_ < t\G()t 2 + t 2 Goti > , 
l ~ ^d^ ' 

The first term Si describes the double scattering with an intermediate NN interaction. It 
is essentially an "optical model " double scattering term. The second term Qi includes 
excitations of nucleons to free continuum states. The expression in Eq.(|7|) is analogous to 
the one in Eq. ([!]), since it contains terms linear in < t > in the denominator. However, 
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it also now contains the prescription for how to calculate the effective distance Rd and 
additionally contains sizable effects from the excitations. Again, the expansion parameter is 
< t > /Rd and it does not need to be small to guarantee the success of Eq.(^), which is also 
applicable if the multiple scattering series is divergent. Higher orders of < t > /Rd in the 
denominator of Eq.(|7|) are obtained by comparing higher orders in Eq.(^) with an expansion 
of Eq.(|7]) with respect to Oi and Si. In this way the next order approximation is obtained 
as 

T * = r^k^s; < 10 > 

where 

S 2 = [S 2 - (SO 2 ] - [fi 2 - (^) 2 ] - [A 2 - ni E a ] , S 2 = <T ^y jVAfT ° rf> , (11) 



< TS* > 



A 5 



< T® d GNN(tiGot 2 +i 2 Goti) > 



< T° d > 



(12) 



One may continue this procedure to include higher powers of < t > /Rd in the de- 
nominator, but for the present problem the precision of the order of 1% has been already 
reached. This was shown in calculations of the 7]d scattering length, |23], f32 |, and optical 
model calculations for rjHe systems, H. This fast convergence rate occurs as a result of an 



almost exact cancellation of the S 2 
also in the scattering region. 



and higher order S n terms. This property holds 



B. Scattering amplitudes, off-shell extension, unitarity 

In this subsection some details of the rj deuteron scattering are discussed, the i]N inter- 
action is introduced and explicit expressions for the multiple scattering integrals are given. 
Next, the question of unitarity is discussed for energies below the deuteron disintegration 
threshold. 

The scattering matrices T n ,i(E),t ri N(E) used for the series summation have the natural 
dimensions of length 2 . To discuss the cross sections it is more convenient to use the scattering 
amplitudes with the dimension of a length. These are denoted by A v d(E), a r?A r(i?), where the 
convention Im A > is used and the threshold values A v d, a V N are the scattering lengths. 
The relations between scattering amplitudes and scattering matrices are 

A vd {E) = -{27i) 2 fi v dT v d{E), (13) 



a vN {E) = -(2tt) 2 /vvV(£)- ( 14 ) 

For the multiple scattering expansion, these relations introduce factors n v d/Hr)N which reflect 
kinematical differences in the propagation of these two systems. 
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Now, some definitions used in this calculation are written down. The integrations in 



Eqs.(lq)-(|12j) are performed in momentum space. To describe wave functions the Jacobi 
coordinates q^N (the relative NN momentum) and p n ( the rj — NN relative momentum) 
are used. However, to describe meson-nucleon interactions within the 3-body system one 
uses other pairs denoted by {%n-iPn)- A useful relation is (Jnn = pn —p^/^- It allows one 
to express the integration volumes dp^dq^N , dpudq^N also in terms of dpNdp^. The energy 
E of the rjNN system is related to the on-shell momentum p v by E = pf ] /2fi v d + Ed where 
Ed is the deuteron energy of — 2.2MeV. For the free propagator one has now 

G = [E NN (q NN ) + E v (p n ) - E]^ 1 5(p n - p v )5(q NN - q NN ). (15) 

The meson-nucleon scattering matrices are t v x[q V N , q' N , E — E(p N )]5(piy — fi N ), where 
the delta-function is used to conserve the spectator nucleon momentum and E{p N ) accounts 
for the spectator nucleon recoil energy. A typical on-shell a^^E) is given in Fig.l. The 
cusp in the real part at the rjN threshold indicates an attraction that is due partly to 
the iV(1535) resonance and partly to other attractive interactions. It is this cusp and the 
attraction related to it, that is responsible for the attraction in the few-body-77 systems, 
leading ultimately to the formation of virtual or quasibound states. The strength of the 
cusp i.e. the value of scattering length a V N is controversial. The one shown in Fig.l comes 



from the phenomenological approach of Ref . p5[ , but a number of smaller as well as larger 



values have been predicted in the literature. For this particular solution one finds the 
effective range expansion 

l/a rjN (E) = l/a riN + -q% N r vN (16) 

with parameters a V N = 0.75(4) + «0.27(3)/m and r v ^ = —1.50(13) — z0.24(4)/m to cover 
the broad energy range in Fig.l. Such a wide applicability of the effective range expansion 
indicates very short range rjN forces. The attitude of this paper is to keep a V N as a semi-free 
parameter, since it results from an interplay of the well established iV(1535) resonance and 
unknown "background" interactions. On the other hand, the energy dependence of a v ^{E) 
is determined mostly by the iV(1535) which generates a negative effective range. Thus the 
value of r V N given above is believed to be weakly dependent on the particular choice of a V N- 
Notice that a V N(E) is small at energies far below the threshold. This allows a perturbative 
approach in this region. 

Although (ijjN is "known" only on the energy-shell one also needs an off-shell extension. 
A separable one is used here: a vN = v ri (q V N)a v jsr(E)v TI (q: ' N ) with a Yamaguchi form factor 



v 



v = (1+g^v/fi^) -1 . The inverse range parameter k v is apparently very large for two reasons. 
Firstly, there is no long ranged meson exchange within the 7]N forces and, as said before, 
the effective range expansion converges rapidly. Secondly, the form factors in any S-channel 
resonance must correspond to the small sizes of these objects. It has to be stressed that 
in the r\d system at threshold one is close to a zero energy binding situation. The results 
are sensitive to every parameter involved including the uncertain k v . A tentative value 
k v = 3.3/m _1 , which corresponds to a natural radius of iV(1535), is used. The sensitivity 
to k v , if it arises, will be indicated. On the other hand, to simplify notation all formulas 
are given in the zero range r\ N force limit. However, the actual calculations are performed 
without this restriction. 
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Within the three-body system, the r/N scattering matrix is to be averaged over some 
energy region, generated by the recoil of the spectator nucleon. So, the rjd scattering length 
given by the impulse approximation formula (R) and normalisations in Eqs. (IT3JT3) becomes 



2 

|0 ( m jp „ \ _ o AV ^ f j^i (=* - m\„ I t? P 



A» d ( Pi ,E jPf )=2^- < / dfr d (p-? i /2)a vN (E--^-)<i> d tf-f f /2) >, (17) 

where pi and p~f are the initial and final meson momenta. Eq.(|TT|) has been obtained with 
initial and final plane waves. In order to extract the S wave contribution an average over 
directions of these momenta is to be performed and this is indicated by the brackets <>. For 
on-shell scattering conditions one has Pi = Pf = p v - This formula is more general, however, 
and provides an off-shell extrapolation for A® d (pi, E,pf) to situations where these quantities 
are not related. The integration in Eq.([L7]) averages the ci v n(E) amplitude over some region 
of recoil energies p 2 /2^m,tiN ■ The extent of this region is determined by momenta involved 
in the deuteron wave function. At the threshold it starts with a negative Ed and extends 
towards negative energies down to about —20 MeV. For positive energies, E covers a range of 
about 20 MeV that includes the r/N threshold cusp. For a limited region of E it makes sense 
to introduce an effective amplitude a vN (E) that has been averaged over these 20 MeV or so 
of the recoil energies. We use Eq.(|T7D to define such an average a v ^(E) and this value will 
be used to calculate higher order scattering terms. In this way the impulse approximation 
amplitude becomes 

A^E^f) = 2^a vN (E) < f dpif> d (p-pi/2)<f> d (p-p f /2) > . (18) 

It is expressed in terms of a deuteron formfactor 

F d {PuP f ) =< J dp4 >d (p-p l /2)4 >d (p-p-* f /2) >, (19) 
which may be also presented in terms of space coordinates as 

F d (Pi,Pf) =< J drexp[ir(pi/2-p f /2))(f) 2 d (r) >= J df j (rp i /2)(p 2 d (r)j ( y rp f /2). (20) 
When the multiple scattering series is summed in Eq.fllDD one obtains the result 

Ar,d(pi,E,p f ) = n s • (21) 

The effect of free continuum r]NN states is, to the lowest order in d V N, given by As the 
integrals in Eq. (|8|) are lengthy we reproduce this formula only in the limit of zero range r)N 
forces 

f) F d (p uPf )J(2n)^ vN [E NN (^) + E v (q + cf)-E] 1 ' 

Now, this equation is used to define an inverse radius of propagation between two successive 
collisions 
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Sl x (p u E,p f ) = a vN ^ <->3B- (23) 

Inverse radii analogous to < - >3B appear in many meson-deuteron models, in particular in 
pion-deuteron calculations. Later we use these to compare approximations and fix potential 
parameters. 

To proceed any further, the nucleon-nucleon scattering matrix tjsrN(q,q', E) is needed. 
It may be expressed by the NN scattering amplitude via a relation analogous to Eq.fllll). 
To generate this matrix a separable Yamaguchi potential Vjvjv — vnn(<i)^nnVnn(<i') is used 
with form factors v^n = n% N /(q 2 + k%n)- The two free parameters k/vat and \nn are 
fixed to reproduce the deuteron binding and one of the following two parameters: the NN 
scattering length or the inverse deuteron radius < - >. The details and motivation for the 
latter choice are presented later. The t^^ enters multiple scattering series of the S n and A n 
type. As the relevant integrals are lengthy we reproduce only the leading term in the limit 
of zero range rj N forces and energy averaged a n N- With this simplification 

2a V N r dp p 2 



Xi( Pi ,E,p f ) = avN [ —^—t NN (E - ) < VfaftVftff) >> (24) 

Fd{Pi,Pf) J {27r) 2 fi 1lN 2fj, NtVN 

where 



The angular average over Pi,p"f must be done to extract the S-wave contribution. Eq.(f24|) 
defines another inverse radius of propagation. 

V x {p i: E )Pf ) = 2a vN -^ < - > 3B • (26) 

The inverse radius <C - ^> differs from < - > defined in Eq. (|23|) above. Approximately, the 
former is quadratic and the latter is linear in the deuteron density and this is indicated by 
the notation. 

The unitarity condition for the on-shell scattering amplitude follows from Eq. (|21|) . Let 
us concentrate on the energy region between the threshold and the deuteron breakup. At 
low energies t NN (E) is dominated by the deuteron pole and it may be presented as 

t NN [E - fi„(p)] Ed + _ E , (27) 

where Nd is the deuteron wave function normalisation and Ed is the deuteron binding energy. 
The absorptive part of t n n comes from this singularity. Below the breakup it guarantees 
the unitarity for the rjd scattering amplitude. To see that, one extracts the absorptive part 
oHnn(E) equal to inN^5[E — E v (p) —E d ] and performs the integration in Eq.(|24|) to obtain 

IrrJ± = 2^F d {p n ,p v )p v . (28) 

The same procedure shows exact cancellation of the absorptive deuteron pole contributions 
to the triple scattering term S2 in Eq. ([T0|) . In addition, an inspection of Eqs.©, ( [P2|) shows 
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that there is no pole in the integrands of the Q n terms. In particular the singularity that 
enters Q\ is related to the cut in Gq. This generates a branch point at the deuteron breakup 
threshold and absorptive contributions at higher energies. Altogether, one can present the 
meson-deuteron scattering amplitude in the form 

as required by unitarity. Below the deuteron breakup threshold and for closed nN channels 
the K v d(E) would be a real analytic function of energy. Above this threshold, the r]NN 
continuum induces an imaginary part to K v d(E). It comes mostly from Qi, much less from 
Si and effects of higher orders are very small. 

To describe final states in the rjd system a half-off- shell scattering amplitude is 
needed. It is obtained by relating pj and E, in the single and multiple scattering terms 
Fd(pi,pf(E)),T,i(pi, E,pf(E)),D,i(pi, E,pf(E)) etc. In the region of our main interest: 
E ~ Ed and pi in the up to \frn~ 1 range one finds a moderate (up to 25%) fall of 
Fdipi, 0). On the other hand the other quantities £i(pi, 0, Ed), 0, Ed) are almost con- 

stant. This occurs since these functions are defined by the ratios of two functions that fall 
down with increasing pi in the same moderate way. The consequence is that the half-off 
shell amplitude may be presented as I square some round brackets in next eq. 

A, nd [ Pl ,E,p f (E)] = ^^W-n^-SW m 

For further use in the description of the rjd system in the process of rjd formation, Eq. (p0[) 
is extended into an operator form 

A°d 

Avd = 1 — ^i(E) — Qi(E) — S 2 (E) ' (31) 

The process of formation and final state interactions requires more elaborate calculations. 
Approximations are done, and to study the value of these we first discuss some standard 
approximations used in elastic meson-deuteron scattering. 



C. The static nucleon and deuteron dominance approximations 

In the literature these approximations appear in many forms. The common element is 
that intermediate states of the nucleon-nucleon pair are reduced to the state of the deuteron. 
Another, related approximation neglects the recoil energy of the spectator nucleon. Although 
these approximations are rather crude this method is useful to discuss limiting situations. 
Below, it is used to indicate a simple interpretation of the terms in our multiple scattering 
expansion. 

The first double scattering contribution is the fli term of Eq. (|22|) . If the nucleon recoil 
Enn is neglected, one obtains at E — 

fifrffc = a nN ^L < - >, (32) 
p vN r 
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where the relevant "inverse deuteron radius" is 



< l >= t df ML£. (33 ) 



r 



The < - > is a standard parameter in low energy pion deuteron scattering [|3~5|1 . It describes 



the meson propagator at zero energy, averaged over the deuteron density. One finds a 
weak but significant dependence of this parameter on the NN potential used to calculate it. 
In particular < = >= 0.448//m for the Paris NN potential [|34j| , f35| and 0.463//m for the 



Bonn potential [[JIJ, P6 |. In this calculation the separable Yamaguchi potential is used which 
generates the Hulthen deuteron wave function. Two sets of potential parameters are tried. 
The standard choice for the potential strength and range reproduces the deutron binding 
and pn scattering length (Yamaguchi I). These are obtained with k^n = 1.41/m _1 and 
= — 0.52/m. However, this choice generates < £ >= 0.559/m _1 , which is too large 
when compared to the superior calculations from the Bonn and Paris potentials. To improve 
the separable model its parameters are fixed to reproduce the deuteron binding and the 
average of < - > obtained with these two local potentials. This requires knn = 0.91/m _1 
and 1/Xnn — — 0.289/m (Yamaguchi II). The np scattering length becomes 6.02/m instead 
of 5.40/m but, for this calculation, the deuteron radius is more important than the np 
scattering length. The results are given in Table 1, where the comparison with the exact 
three-body calculation is made. A factor of two difference is found. It comes from the large 
nucleon recoil entering the propagator in Eq.(|22"D and the deuteron binding. The latter is 
introduced into the three body calculation by the condition E = E d ( not E — 0) at the 
threshold. 

Now, we proceed to study the other double scattering term Ex. At very low energies the 
deuteron dominance is expressed via 

GNN{q > E > q) ~ E- Ed -E v ( P y (34) 

This formula greatly simplifies all the multiple scattering integrals. To present the argument, 
the threshold scattering (E = Ed) is considered. The dominant term Si becomes 

£f atic = 2a vN ^- « - », (35) 

where another inverse deuteron radius is introduced 

« I >>= f m p Mr)W)\ (36) 
r J \f/2-r'/2\ 



in close analogy to the more general three-body definition of Eq. (pb|) . The factor of 2 in 
the denominator is essential, since Ylf aUc corresponds to the double scattering amplitude 
generated by an optical potential 

V stauc {r) = _^[ 2a ^][2Vg(2r)]. (37) 

The strength of this potential 2d V N is normalized to two nucleons. The profile <f) d (2r) is 
"compressed" as it refers to the center of the deuteron, hence the argument 2r and the 
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normalization factor 2 3 . y statlc gives the correct expression for the single scattering am- 
plitude but misses badly the double scattering term. In the static approximation one ob- 
tains « 1/r ^>= 0.83/m" 1 while the correct three-body value is <C 1/r ^>3#= 0.32/m _1 
(Hulthen-Yamaguchi I wave function and zero range r]N force, see Table 1). This large 
difference is, again, related to the nucleon recoil in the three-body propagators. Now the 
propagator enters Eq . (p4j) twice and this makes the static approximation even worse than 
in the < 1/r > case. 

For scattering from a system of A particles, a "remedy" for this inconsistency has been 
found long ago by substituting the AA factor by A(A — 1). For the deuteron this leads to 
a reduction of the double scattering term by one half. One can see in Table 1 that such 
a procedure brings the static optical potential model into reasonable consistency with the 
sum of Si and Qi terms. 



D. The pionic channels 

So far the effect of pionic channels has been hidden in the absorptive part of the elastic 
i]N scattering amplitudes. These describe processes related to a single nucleon. In addition 
there exist processes when an intermediate tt meson is exchanged between two nucleons. 
The effect of such intermediate states is expected to be small for two reasons. First, the 
transition from the rj to the 7r meson involves a 300 MeV energy release. This energy is shared 
between the three participants N, N and ir. The bulk of the NN interactions occur at fairly 
high energies where the NN scattering amplitude is small. Second, at these energies the 
intermediate state propagator oscillates in space at distances much shorter than the deuteron 
radius. This suppresses the multiple scattering integrals. However, since the NNrj system 
is close to binding even small effects may matter. On general grounds one could expect that 
the NNn channel which has a lower threshold would contribute some repulsion into the 
NNrj channel which has a higher threshold. That could reduce the chances for binding in 
the NNrj system. 

With the explicit pions the series for the 7/-deuteron scattering matrix T n d is supplemented 
by terms 

AT vd = hiG h' 2 + h 2 G h'i + (fn + h 2 )GW(A'i + h' 2 ) + ... , (38) 

where the scattering matrices h describe the rjN — > nN transitions and the h' describe 
ttN — > rjN transitions. These are related to the scattering amplitudes via Eq . fll4|) . However, 
the kinematic factor is now y^jv/^riv. The small pion mass reduces the effect of these factors 
as the relevant ratio fi n N/^nd is close to unity. 
The pion exchange scattering is described by 

^(Pi,E,p f ) = _ °l N f N M m*), (39) 

where I(m n ) has the form of the integral in Eq.(^) but with the mass of the rj meson entering 
this integral being changed into the mass of the ir meson. The expression Eq. ( p9| ) corresponds 
to the free propagator part of the series in Eq.(^). The other terms are described by a 
formula similar to Eq.(|24) 
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On 2 

ZJ(p h E,p f ) = _ -Sf* /Km,), (40) 



where I(m n , m^) is given by the integral in Eq.(|24T) but with the mass of the n meson being 
replaced by the mass of the ir meson. 

A few numbers calculated at the rjd threshold are now presented. These correspond to 
the older version of the rjN model of Ref. p5| . The scattering parameters are : = 
0.57 + z0.14/m and cl^^n — 0.21 + iO.lbfm, while the double scattering integrals are 
I(m*) = (-3.5 + il4.9)10- 2 /^ -1 and I(m*,m*) = (6.89 - i3.93) lO" 3 ^- 1 . Thus the Ef 
is a negligible 1/pm correction to the multiple scattering. The other term is significant. 



Corresponding numbers for the elastic rjN channel may be found in Ref. |23| and in Table I. 

The calculation of without the ir exchange effect and the precision obtained with 
the multiple scattering theory (MST) can be estimated by comparing with the solution of 
the Faddeev equations (F) obtained in Ref. |32|] with the same input. For small a V N these 



methods give essentially the same answer. At larger values, close to the binding situation 
differences arise. These are given in Table II. 

The MST method slightly underestimates the real parts and overestimates imaginary 
parts. 

The dominant source of these differences is related to the use of energy averaged am- 
plitudes a V N in the second and higher orders scattering terms. No attempt is done here to 
correct for that. Also there are higher order effects and slight differences in the effective 
range corrections involved in these two models. 

Our best result for the A vd scattering length obtained with the elastic r/N amplitude 
of Ref. |p5| and the Yamaguchi II model for the NN interactions is A^ = 2.31 + il.73fm. 



This number should be compared with A vd = 2.61 + il. 72 fm obtained in Ref. [33| from the 



Faddeev equations, using the same a V N but a somewhat superior separable NN potential. 
The results given above do not include the pion exchange effects. If the pion exchange is 
introduced we obtain A v d = 2.27 + il.68fm. The change is rather small and corresponds to 
a repulsive effect. Qualitatively, such a behaviour has been predicted in Ref. |33| but the 



size of the effect obtained there is much larger and is not reproduced here. The difference 
may be related to the small values of cl^^^n inherent to the K-matrix model of Ref. [|25 



Eq. (39) indicates that VL\ is composed of three uncertain complex numbers and may be 



unstable with respect to the input values of a V N,irN an d a V N- Here, the averaged propagator 



/(m,,-) has been calculated with non-relativistic kinematics as used in Ref. [p3f1 . A relativistic 
expression for the pion energy would make the pion exchange effect even smaller. 

Numerical results for the rjd scattering amplitude A^E) and the rjd elastic cross section 
are given in Figs. 2a ,2b. The impulse approximation is valid at energies E > 15MeV above 
the threshold. Multiple scattering effects enter at lower energies : in the intermediate region 
of E fa lOMeV it is Q\ that dominates while for lower energies the main contributor is Ej. 
The real three-body exotics is seen only very close to the threshold for E < lMeV + E d . 
The effect is strong and very narrow. At present a direct way to see it experimentally is not 
available, and one has to study inelastic processes that involve the nd system in final states. 
One reaction of interest is pn — > dn, with a related one being the pp — > ppn. In the next 
section the deuteron reaction is discussed. 
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III. FINAL STATE INTERACTIONS IN THE i]D SYSTEM 



The amplitude squared for the pn — > rjd process extracted from the CELSIUS experiment 
is given in Fig. |j. Three regions may be specified : 

(1) The "asymptotic" region of E > 20MeV, where the final state interactions are small, 
can be described by the impulse approximation. 

(2) The "interference" region of E m lOMeV, where the experimental amplitude indicates 
a shallow minimum. 

(3) The region of a "virtual three body state" located below the deuteron breakup thresh- 
old. 

In this section we find F v d, the amplitude for the pn —>■ drj reaction. It consists of several 
terms that differ in their physical content. In order to link this amplitude to the standard 
description of two-body final state interactions, it is presented in the form 



A vd {E) ] 
Rd(E) 



^ d = A c [l + ^^}. (41) 



This expression consists of the on-shell 7/-deuteron scattering amplitude A v d(E) and a radius 
Rd(E). For simple models, like the static nucleon approximation, a simple interpretation 
and a closed formula for l/R^E) are found. The radius comes out as an interplay of other 
radii: the n meson source radius, the nN force range and the deuteron radius. For scattering 
on static as well as interactive NN pairs this radius makes sense only at very low energies. At 
energies exceeding about 10 MeV it becomes strongly energy dependent. The presentation 
of FSI in terms of Eq . ([41]) is still possible but the simple physical interpretation is lost. 

The three regions indicated in Fig. |] are specified for three different purposes. The main 
emphasis is put on the low energy enhancement and properties of the three-body state. 
There, it is the ratio A v d/Rd that represents the strength of the enhancement. In order to 
check on A v d one needs full control of 1/Rd- The interference region indicates an energy 
dependent phase of l/Rd(E) that together with the phase of A vd (E) produce a minimum. 
Finally the cross section at higher energies is used to fix A c which is an amplitude for 
coalescent formation of the deuteron. To calculate it, one needs a specific model of the 
meson formation. This question is discussed in the first subsection. 



A. Meson formation amplitude 

The formation of an n meson in nucleon-nucleon collisions involves large momentum 
transfers and the mechanism of this process is uncertain. It is believed to be dominated by 
an intermediate p meson exchange, although other mesons may also contribute |5H, jj, jj, 



4"2|| , |f43| . In particular the recent calculation, ||44j| , stresses the role of n meson exchange. 



All models assume three stages for this process, namely, NN — > iViV(1535) — > NNij, where 
the intermediate stage is dominated by the iV(1535) resonance. The first transition involves 
large momentum transfers, whereas the second is a low energy process. 

This paper aims at the final state effects and the formation amplitude is described in a 
phenomenological way by a function of initial and final momenta A^ orm . What is essential 
for our study is to find a functional form for this amplitude and pinpoint the effects (if any) 
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that may generate a rapid energy and momentum dependence. Such a dependence would 
be equivalent to a long range space structure of the rj meson source and this structure is to 
be viewed on the scale of the deuteron radius. 

The 77-meson production model is illustrated in Fig. [3[ It consists of a short ranged 
meson exchange part affected by initial NN interactions. As compared to the deuteron 
radius it is obviously a short range process, since it requires large initial NN cm. momenta 
Qnn ~ 3.5/m" 1 and comparable momentum transfers. This part of the production process 
may be described by an amplitude (q 2 + /i 2 ) -1 , where q t is the momentum transfer and fi is 
the mass of the exchanged meson. Now, the momentum transfer may be expressed in terms 
of the initial and final momenta q t = —Qnn — Pin, where the final nucleon momentum is 
Pin = Qnn — Pri/2. For the final states of interest, Qnn 3> Pin an d this amplitude is almost 
constant. Projected into an S wave it generates a momentum dependence (Qnn+I^+Pin) 1 - 
As might have been guessed already from the uncertainty principle, the range involved 



l/yQ 2 NN + n 2 is less than 0.35/m even for 7r-meson exchange. 

However, there is another part of the formation mechanism, related to the propagation 
of the iV(1535) resonance relative to the spectator nucleon. The propagation range is not 
negligible on the scale of the deuteron radius. It may be described by a function 
G r = [A — iT(E)/2 + E re( \~ l , where A = E r — Mn — M rj — E, the nucleon recoil energy 
Erec — Pin /2aVv,tv; an d E r — iT(E)/2 are the resonance energy and width. It is the recoil 
energy that introduces a sizable momentum dependence in this propagator and generates 
a long range structure of the formation process. Altogether, the formation amplitude is 
assumed to be of the form 

/1/27T 2 

«; + (pin) 2 

where k s is the inverse radius of the rj formation source. The Fourier transform of this 
amplitude with respect to the nucleon momentum p iN becomes 

Af 

A form (r) = —exp(-K s r). (43) 
r 

In general k s , given by the relation n 2 = 2 / u 7? 7v,A r [A — iT(E)/2], is energy dependent but our 
main interest is the rjN threshold. The resonance position and half width of 1535 — 2148/2 
MeV are taken from Ref. ||25|| . At the threshold, the decay of the iV(1535) into the r]N 
channel becomes negligible and the actual width there is about 60MeV. These parameters 
yield a complex value of n s = 1.40 — z0.39/m _1 . Strength of the meson formation amplitude 
is a free parameter. It is used to fix the calculated pn — > ijd cross section to the 
experimental one in the "asymptotic" region. 

B. Final state wave function 
The full amplitude F v d for reaction pn drj is given by an integral 

F vd = J dqdpAf° rm (q- p/2)^- d *(q,p), (44) 



15 



where the nucleoli momentum is expressed in terms of Jacobi NN and r) momenta. In this 
equation is a wave function for the final three-body system that fulfills the ingoing wave 
boundary condition. To find it, let us construct the multiple scattering series for the wave 
function similar to the series for the scattering matrix. From the Faddeev equations one 
obtains 

<S> vd = [1 + G f vd + G NN f vd }(f) d ^ v , (45) 

where T nd is the scattering matrix operator given by the series in Eq. [|. Now the reaction 
amplitude contains three distinctly different terms 

F vd = A C + A ex + A vs (46) 

that correspond to the three terms in the square bracket of Eq. (|4"5|) . The first one is a direct, 
coalescence transition from the source to the deuteron : 

A c = J dqdpAf° rm (q- p/2)MQ)MP)- ( 4 ?) 

For an S-wave free meson final state of momentum p v and the specific choice of the formation 
amplitude in Eqs. fl4"2|j43|) the coalescent amplitude is given by the integral 

A c = J dfAf° rm (r)Mr)Urp ri /2). (48) 

The formation amplitude is short ranged in comparison to the deuteron radius. Thus the 
coalescent formation amplitude involves <fi d (r m 0) and in the region of interest it depends 
weakly on the meson momentum. 

The second term in Eq.(f46|) describes a part of the rj meson final state interactions 

A ex = J dqdpA^G^p^M,. (49) 

The superfix ex indicates that this term contains an rj meson exchange diagram. Here the 
rjNN system propagates freely from the formation source and this process is described by 
Gq. Next, the scattering series begins with the meson-spectator- nucleon collision. Thus the 
meson produced on one nucleon is exchanged and interacts with the other nucleon. After 
that the deuteron is built either in the coalescent manner or in subsequent NN interactions. 

The last term (A vs ) in Eq . (f46"|) is our main interest. It describes the final state scattering 
process that begins with the NN interactions. These are the strongest interactions in the 
rjNN system. For low energies the corresponding propagator Gnn is dominated by the state 
of the deuteron. The related amplitude is 

A vs = J dqdpA^ m G NN {q,p)f nd M v . (50) 

On top of the deuteron described by G^n the system builts up the rjNN virtual state (vs) 
and this process is described by T vd . It is to be seen as a narrow enhancement on the 
background of the standard coalescence term A c . 

Before discussing the r]NN collective state let us return to the meson exchange process 
and the related amplitude A ex of Eq. (f49"D . The final state scattering begins with the rj- 
meson propagating from the source to the other nucleon and in leading order it is described 
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by Got V N- To this order one can rearrange the scattering series in Eq. (^) and the final state 
wave function to the form 

$ vd = (1 + G i vn )[l + GWFi J d]0dV'»j> ( 51 ) 

which reduces the FSI to the coalescent (A c ) and virtual state (A vs ) amplitudes. However, 
now the the r\ formation amplitude is redefined to A^ orm [l + G f o^Jv]- This formation process 
is of an extended range in space but is still much shorter than the deuteron radius. The 
initial meson exchange process involves nucleons correlated at short distances and thus 
requires large momenta and a large recoil. The relevant energies in the rjN system extend 
far below the threshold. There the a V N amplitude is not known. According to the model 
shown in Fig.l it is expected to be real and small. Following this we terminate the exchange 
process at this level. The result of this analysis is that the unknown strength of the i] source 
is renormalized and its size may be somewhat larger than that predicted by the iV(1535) 
resonance model. In the numerical calculations a plausible value of k s = 1.5/m -1 is used to 
describe the three factors that contribute to the size of the i] meson formation source. 

The task now is to find the coupling of the virtual state to the final i] - deuteron channel, 
and again the method of Sect. II is used for this purpose. We begin with the simplest 
calculation. 



C. The deuteron dominance approximation 

This approximation has been systematically used in all calculations of the rj final state 
interactions on light nuclei. In this section we discuss its applicability in the deuteron case. 
Ultimately it is shown that the use of static nucleons is a bad approximation. On the other 
hand, this method is useful to discuss some approximations. We take advantage of that. 

At very low energies the deuteron dominance of G^n is expressed by Eq. (|34j) . This 
formula inserted into Eq.(^U[) gives A vs in terms the scattering matrix T vd . One finds 

A- = J dqdp Af™(q- p/2)Ml) E^E^-E - (52) 

The right hand side of this equation expresses A vs in terms of the half-off-shell T vd (q, E,p v ) 
scattering matrix. The problem that arises at this stage is the off-shell continuation of 
T vd . As discussed in the previous section the main dependence of the off-shell extrapolation 
comes essentially from the numerator of Eq.(^) i.e. the deuteron formfactor of Eq. (|19|). 
Within this approximation it is easy to extend the scattering to the case of intermediate 
mesons which are described by plane waves. The virtual-state amplitude in Eq.( may 
be presented in the coordinate representation. With the use of the intermediate meson 
propagator 

exp(ipR) exp{ip v R) 

P {2^f Nd [E d + E v {p)-E] R 1 } 



and Eqs. (p~8| , p~9|) one obtains 
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A- = j d p d? A f-(r')Mr') ^IrZi r '/ 2 ^(r) j (r Pt? /2) . (54) 

Fd{Prj,P n )J \r/2-r'/2\ 

The basic amplitudes A c and A 1141 are now easy to calculate in two extreme cases : of a zero 
radius and a large radius meson source. For the zero radius source given by A^ orm {r') = 
A* 5{r') one obtains the coalescence amplitude 

A c = At Mr = 0), (55) 

and the total amplitude which includes final state interactions 

F« = A° + A" = A°[l + ^]&), (56) 

where the inverse interaction radius is defined as 

drMr) jj- jo{rp v /2). (57) 



R d (E) F d (p v , Pr] )J rdW r/2 

With the zero range source and static nucleons the radius Rd{E) is due entirely to the 
deuteron size. For energies close to the threshold Eq.(|5"TD relates it the previously defined 
static inverse radius 

-L = 2<->+i P „ + 0(p 2 7 ). (58) 

Numerically one obtains Re I/R4 = 1.12 fm~ l for the Hulthen I deuteron wave function. 
The expansion in Eq.(|58|) adds a simple imaginary part to < - >. For static nucleons this is 
an exact expression as may be easily checked with Eq. (|20D for the deuteron formfactor. For 
the FSI in two-body systems described by a potential, such an imaginary term is due to the 
imaginary part of the outgoing wave exp(ipr)/r. It gives rise to an interference effect in the 
square-bracketed term of Eq. floTf) . By splitting the FSI amplitude into real and imaginary 
parts one obtains 



* rjd 



A c I 2 



1 1 

[1 + ReA r]d < - > -p^ImA^} 2 + [ImA^ < - > +p ri ReA vd } 2 



(59) 



As the amplitudes in this equation are positive, one notices that, with increasing momentum 
p v and constant amplitudes, the first term in Eq. ( [59"D decreases and the second term increases. 
This leads to a minimum in the amplitude. Such a minimum is indicated by the experimental 
data in Fig. |] in the region above the enhancement, where the meson deuteron scattering 
amplitude becomes more stable. 

The static nucleon approximation is convenient for the other simple limit of a large 
meson source. If A^ orm (r') = A^0rf(r), i.e. the source is as large as the deuteron, then the 
low energy expansion is given by formula (j36D 



=< _ >+ ip +0(pj). (60) 
Rd r 1 

The radius R d given by this equation is sizably larger as one finds jj- = 0.83/m -1 . The 
two examples given here indicate the sensitivity of the final state interactions to the source 
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radius. For a more realistic source determined by the propagation of the N* resonance via 
Eq.© we find Re±- = 0.90/m" 1 , at the threshold. 

However, static nucleons are a rather poor approximation. The results are given in Table 
III, but as may be guessed from Table I the values of l/Rd are likely to be overestimated by 
a factor of 2 to 3. To improve this situation it is necessary to analyze the dynamic deuteron 
formation via the NN interactions. This is studied in the next section. 



D. A dynamic formation of the deuteron 

The main approximation behind Eq.(^) is a reduction of three body propagator to a 
quasi two body propagator. On the other hand, a correct expression for the AW propagator, 
that follows from the definition in Eq.(|3]) is 

n ( i T7i\ v NN (q)t NN [E - E ri (p)]v NN (q') 

GNN{q ' q ' E) ~ [E NN { q ) + E ri { P )-E][E NN { q >) + E ri { P )-EX (61) 

For the threshold energy E = Ed the denominators in the three body propagator involved 
in this equation have no singularities. Together with the form factors vnn these generate 
the deuteron wave function 

, f v N d v NN (q) 

Mq) = W^j - EA <62) 

provided the energies E v (p) in the denominators are dropped. The normalization factor 
results from the residue oft NN at the deuteron pole. Thus the static nucleon approximation 
of Eq.(^l]) consists in dropping the E v (p) in the propagator and reducing the t NN to a 



bare singularity of Eq . (|27|) . These two approximations, in part, tend to compensate each 
other. Roughly, the two energies in the denominators E n (p) and Enn(q) are comparable two 
energies in the denominators E v (p) and Enn{<i) are comparable and the deuteron dominance 
overestimates Gq by a factor of 2 2 . This happens because the deuteron is a weakly bound 
object and the binding energies in the denominators are small in comparison to the recoil 
energies. On the other hand by limiting tjvjv to the pole term one neglects a sizable attraction 
that exists in the AW system at large negative energies. The pure pole term of Eq.(P?D 
underestimates t/viv by a factor of 1/2. The net effect for G^n is a reduction to about a 
half of its static nucleon values. The details depend on the intermediate momenta involved 
in this process. 

In the dynamic approach the virtual state amplitude A vs is given by an up-dated version 
of Eq. (|52|) , where now the deuteron pole expression for Gnn given by Eq.(|3~4"D is substi- 
tuted by the full expression of Eq.(|6l|). In addition, the half-off-shell r\ deuteron scattering 
amplitude given by Eq.(|30|) is used. The full A vs is given by an integral similar to the in- 
tegral in Eq.(|24D for the Si. However, in Eqs.([24]),(^) one of the deuteron wave functions 
(f)d(q — should be replaced by the source function Af° rm (q — ^p). In the limit of the 
deutron domination this procedure returns back to the simple equation ([52]) . If one wishes 
to present A vs in the form of Eq.([TI|), then the inverse radius becomes a ratio of two 



integrals. At the threshold, the effects of the three body propagator and source size reduce 
this quantity to R ^ = 0.46. This number is more than a factor of 2 smaller than that 
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given in the zero range source, deuteron dominance model. Some other values calculated 
for two NN potentials and two source ranges are given in Table III. It is also found that the 
energy dependence of Rd{E) is more moderate than that given by Eq.([|7]). This reflects the 
significance of the attraction in t^N at negative energies below the deuteron pole. 



IV. RESULTS 

The elastic r]d cross section is plotted in Fig. 2a for two values of the inverse scattering 
parameters < - >. These parameters differ by 15% only, but the effect on the cross section is 
significant. It reflects proximity of the r\d quasibound state. A similar sensitivity is displayed 
by the cross section on the value of the 7]N scattering lengths. A few possibilities are plotted 
in Fig. 2b. These differ by the input scattering lengths a V N and effective ranges. Since the 
cross section is not measured, it cannot be a test for a V N- The point of interest is the relation 
of the peak in the cross section to the position of the quasibound singularity in the complex 
p v plane. This singularity is found as the position of zero in the denominator of Eq . (pip . 
An approximate solution for a V N = 0.75 + i0.27fm case is p vs = —0.09 — iO.Odfm^ 1 and it 
is located in the third quadrant of the complex plane. This corresponds to a virtual-state, 
that is analogous to the spin singlet proton-neutron state. Coupling to the pionic channels 
pushes the singularity from the negative imaginary semi-axis into the third quadrant of 
the p v plane. This singularity is far from the region where quasi-bound r]d states might 
exist. Those would be analogous to the deuteron. To enter the quasi-bound state region 
with the model discussed here, in particular with the subthreshold amplitudes shown in 
Fig.l, one would need Re a V N > 1.2 fm. The K- matrix model based on two body data 
allows scattering lengths larger than 0.75/m, up to about 1.05/m. The corresponding 77 — d 
amplitudes are given in Fig. 2b. However, as discussed below the large values close to lfm 
seem to be excluded by the rjd final state interactions. 

In Fig. f| the CELSIUS data [0 are plotted and compared with calculations. The cross 



section around the 10 MeV excess energy region is used to fix our normalisation constant 
A? . The choice of fixing point is not particularly relevant, provided it is done above the 
low few-MeV energy region. However, as is seen from the figure, our calculations are on 
the lower side of the experiment for E > 15MeV and on the upper side for E < 3MeV. 
The 9 MeV point chosen helps to make a good overall fit to the data. At low energies the 
experimental energy resolution is vital. The experimental points refer to the middle of the 
corresponding energy bins. No attempt is made to fold the calculated cross section over the 
experimental energy resolution. At low energies and in the enhancement region this cross 
section is a check on the r)N model. The data and the results calculated with the Yamaguchi 



II NN potential and r]N interactions of Ref. [25] are consistent. The same data with the 



experimental errors enlarged by the beam energy uncertainty (0.3 - 0.5 MeV) could also 
accomodate the results obtained with the Yamaguchi I NN potential. Large values of a V N 
would not be allowed, at least with the subthreshold behavior generated by this K-matrix 
model. 

The position of the minimum in the scaled cross section comes out as an interference 
effect. It depends on details of the energy dependent rjd scattering amplitude. However, 
with the present precision of the data and uncertainties of the input, it does not present a 
test for the consistency of the data and calculations. 
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TABLES 



TABLE I. Inverse deuteron radii ( propagator at zero energy averaged over deuteron wave 
functions) calculated with different deuteron wave functions and/or different propagators for the 
rjNN system. The <> entries correspond to free continuum states: <> static nucleons, <>3B 
3-body propagator. The <CS> entries correspond to intermediate NN interactions : <C3> projection 
onto the deuteron state, <C3>3_b complete NN interactions. The three-body entries do not include 
effects of a finite rjN force range. With k v n = 3.3/m -1 these should be reduced: <>3_b by 18 % 
and ^C3>3.b by 13%. All entries are in fm^ 1 

Paris Bonn Hulthen II Hulthen I 

< i > 0.449 0.463 0.456 0.559 

< p >sb 0.197 0.254 
«i» 0.721 0.831 

< i >3B 0.280 0.324 



TABLE II. Comparison between the Multiple Scattering Theory (MST) used in this paper and 
the Faddeev(F) results for — see Ref. pal 



a^(fm) yl^ (fin) (MST) A vd (fm)(F) 

0.44+i0.30 1.01+il.50 1.03+il.49 

0.62+i0.30 1.65+i2.41 1.71+i2.34 

0.888+i0.274 2.37+i5.79 2.65+i5.48 



TABLE III. Inverse radii < -i- > used for the description of r\d final state interactions at zero 



energies. All entries are in fm 



Method 7] source Hulthen II Hulthen I 

Deuteron Dominance Zero range 0.91 1.12 

Resonance dominated 0.75 0.90 

Deuteron Formation Zero range 0.504 0.555 

Resonance dominated 0.465 0.528 
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FIGURES 

FIG. 1. The r]N — * r]N amplitude a^jy in fm as a function of the center-of-mass energy Eq.m. 
in MeV. Solid line for Re a^jv and dotted line for Im a v ^. 
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FIG. 2. \A v d(E)\ 2 in fm 2 as a function of E(MeV) the rj—d kinetic energy in the center-of-mass 
system . Figure a) shows the results for Re (2^=0.75 fm but with two different NN potentials : solid 
line - Yamaguchi II and the dashed line Yamaguchi I. Figure b) shows the results for Yamaguchi 
II and Re a V N = 0.75 fm (solid line), 0.87 fm (dotted line), 1.05 fm (dashed line) and 0.21 fm 
(long-dashed line) from |25| |29|]. 
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FIG. 3. The ^-formation reaction 
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FIG. 4. The reduced cross section r = a(pn — > i]d)/p v (E) from Ref. [11 normalised to unity 
at 9 MeV. The theoretical equivalent \Fr, d (E)\ 2 /\F vd (9 MeV)| 2 is calculated for Re a r)N = 0.75 fm 
(solid line), 0.87 fm (dotted line), 1.05 fm (dashed line) and 0.21 fm (long-dashed line) from [25] 
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